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Abstract 

We investigate two-site electronic correlations within generalized Hub- 
bard model, which incorporates the conventional Hubbard model (param- 
eters: t (hopping between nearest neighbours), U (Coulomb repulsion (at- 
traction)) supplemented by the intersite Coulomb interactions (parame- 
ters: J W (parallel spins), (antiparellel spins)) and the hopping of 
the intrasite Cooper pairs (parameter: V). As a first step we find the 
eigenvalues E a and eigenvectors \E a ) of the dimer and we represent each 
partial Hamiltonian E a \E a }(E a \ (a = 1, 2, .., f6) in the second quantiza- 
tion with the use of the Hubbard and spin operators. Each dimer energy 
level possesses its own Hamiltonian describing different two-site interac- 
tions which can be active only in the case when the level will be occupied 
by the electrons. A typical feature is the appearence of two generalized 
t — J interactions ascribed to two different energy levels which do not 
vanish even for U = J^ 1 ' = j' 2 ^ = V = and their coupling constants are 
equal to ±t in this case. In the large U - limit for = j' 2 ' = V = 
there is only one t — J interaction with coupling constant equal to At 2 / \ U\ 
as in the case of a real lattice. The competition between ferromagnetism, 
antiferromagnetism and superconductivity (intrasite and intersite pair- 
ings) is also a typical feature of the model because it persists in the case 
U = J (1) = J (2) = V = and t 0. The same types of the electronic, 
competitive interactions are scattered between different energy levels and 
therefore their thermodynamical activities are dependent on the occu- 
pation of these levels. It qualitatively explains the origin of the phase 
diagram of the model. We consider also a real lattice as a set of in- 
teracting dimers to show that the competition between magnetism and 
superconductivity seems to be universal for fermonic lattice models. 

PACS numbers: 71.I0.-w, 71.10.Ca, 71.10.Fd, 71. 20. Be, 71.27.+a, 
71.28.+d, 73.21 La, 78.67 Bf 
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1 Introduction 



Model Hamiltonians which serve to describe electronic subsystems of solids are 
formulated on the quantum mechanical basis which takes into account interac- 
tion processes between electrons, most of them basing on the Hubbard model 
and its generalizations (cf e.g. Refs [l]-[96]). The ground state properties of 
such models has been investigated in Refs [85]- [92]. Each model Hamiltonian 
of this type desribes, however, as a rule many unknown, competitive, elec- 
tronic correlations (two-site, three-site, etc.) which are normally invisible in 
the original model but the knowledge about their existence is very important 
because just these correlations determine also an area of physics where a given 
model Hamiltonian can really be applied. To find the electronic correlation 
within a given model Hamiltonian H let us assume that we can exactly solve 
the Schrodinger equation H\E a ) = E a \E a ). Having to our disposal the calcu- 
lated eigenvalues E a and eigenvectors \E a ) we can use the equivalent form of the 
given Hamiltonian H = ^2 H a where the set of commuting partial Hamiltonians 

a 

H a = E a \E a }(E a \ we can represent in the second quantization for each energy 
level separately introducing here Hubbard and spin operators. In this way to 
each energy level of the system E a a partial Hamiltonian H a can be ascribed. 
Each part H a contains many important, competitive interactions, active in the 
case when the level will be occupied by the electrons. This simple idea cannot 
be unfortunately applied in a general case because we cannot exactly solve the 
mentioned Schrodinger equation. It is, however, possible to do it exactly and 
analytically in the case of a dimer described by the generalized Hubbard model 
(see Sec. 2) to show that two-site electronic correlations resulting from this 
approach describe the competition between magnetism and superconductivity. 
We show that this competition is an universal feature of all electronic lattice 
models containing hopping term and it takes also place in the case of a real 
lattice. This result seems to be very important also in the case of quantum dots 
and nanostructures (cf e.g. Refs [97]- [100]). 

2 The Model 

The generalized one-band Hubbard model belongs to a class of fermionic lattice 
models widely used in the solid state physics. This model has been primar- 
ily applied to explain magnetic and transport properties of transition metals, 
their compounds and alloys, including also insulator-metal transitions (cf e.g. 
Refs [1]- [16] and original papers cited therein). After further generalizations the 
model has also been applied to fluctuating valence systems and heavy fermions 
(Anderson-like models, Refs [17]-[18] (see also Ref. [13] for a review)), liquide 
3 He (see e.g. Refs [19]-[21]) and fullerenes (cf e.g. Refs [22]-[24]). A special 
attention has been, however, paid in recent decades to the theory of high-T c 
superconductivity (cf Refs [25]- [58]) where the extended Hubbard model (neg- 
ative U model) has widely been used (see e.g. Ref. [33] for a review). Another 
interesting model also used in the context of superconductivity is the KPPK 
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model, formulated in Refs [59], [60] (see also e.g. Refs [61]-[64]). In the present 
paper we consider the extended Hubbard model supplemented by the hopping of 
the intrasite Cooper pairs (KPPK interaction). The Hamiltonian of this model 
has the form 



(1) 



The indices enumerate the lattice points (Ri, Rj), tjj is the hopping inte- 
gral, U denotes the effective intrasite Coloumb interaction, and (gener- 
ally, not necessary equal) describe the effective intersite interactions, all of them 
resulting from the original intrasite and intersite Coulomb repulsion which can 
be modified by polaronic effects (see e.g. Ref. [33] for details) and therefore 
U, j( 1,2 ' can be treated here as positive or negative parameters. The last term 
in (1) is responsible for the transport of the intrasite Cooper pairs (Refs [59], 
[60]) with the coupling constant V. The model (1) cannot be solved exactly in 
a general case. We can, however, consider a special but nontrivial case of two 
interacting ions (a dimcr problem) which posseses exact, analytical solution. 
Thus, let us start with the dimer Hamiltonian, resulting from the expression 
(1). It has the form 

H d = -*E( c t<T C 2,<T + cJ (T ci i(J ) + [/(ni, T ni ii +n2, T n2,j.) 

E n ha n 2 , a + J {2) E n ll<r n2,-a (2) 

a a 

-V(cTjcT tl c 2 ,ic 2 ,i + c+ T c+ iCu ci, T ) 

where t 1>2 = t 2 ,i = -t, J^f> = J^f 1 = J (1 < 2) and V 1<2 = V 2 .i = V. We start 
from the Fock's basis \n\^, 711,4; 712,4, 712,4.) { n i,<y = 0, 1; i = 1, 2; <r =|, J.) and we 
find the exact solution of the dimer eigenvalue problem (Ho\E a ) — E a \E a )): 

E 1 =0; |£i) = |0,0;0,0>, 

E 2 = -t; \E 2 ) = ^(|1, 0; 0, 0} + |0, 0; 1, 0)), 

E 3 =t; |^> = ^(|1,0;0,0>-|0,0;1,0», 

Ei = -t; \E 4 ) = ^(|0, 1; 0, 0) + |0, 0; 0, 1)), 

E 5 =t; |S B > = ^(|0,1;0,0>-|0,0;0,1», 

E 6 = J (2) ; \E 6 ) = ^(|1,0;0,1) + |0,1;1,0)), 

E 7 = U + V; \E 7 ) = -^(|1,1;0,0)-|0,0;1,1)), 

E 8 = C + u-v+jiv . |£7 8 ) = a+ (|i j i ;0 ,0) + |0,0;l,l))-a_(|l,0;0,l>-|0,l;l,0», 

E 9 = -C + u-v+J™ . \e 9 ) =a_(|l,l;0,0> + |0,0;l,l» + a+(|l,0;0,l>-|0,l;l,0», 

^io = J (1) ; £7 10 > = |1,0; 1,0), 

£n = J (1) ; \E n ) = 0, 1; 0, 1), 
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E 12 = t + U + JW + JP); \E 12 ) = ^(|0, 1; 1, 1) + |1, 1; 0, 1)), 

E13 = ~t + U + J« + J( 2 ); |£7 13 > = ^(|0, 1; 1, 1) - |1, 1; 0, 1)), 

E 14 = t + U + J« + J( 2 ); \E U ) = 75(11,0; 1, 1) + |1, 1; 1, 0}), 

E 1B = -t + U+ J« + JW; \E 15 ) = ^(|1,0; 1, 1) - |1, 1; 1, 0)), 



E ie = 2(U + JW + J( 2 )); \E W ) = |1, 1; 1, 1) 

where 



(3) 



C=^(^2>) +4*2, a ± = yi±^^. (4) 
In the following we apply Hubbard and spin operators (cf e.g. Ref. [1]) 

5f = 2Kt- n i,i)' = °.>°*^. 



(6) 



and we use the equivalent expression for the dimer Hamiltonian (2) 

16 

a=l 

where P a — \E a )(E a \. Each product E a P a in the formula (7) where we insert 
E a and \E a ) from the formulae (3) can be rewritten in the second quantization 
with the use of the Hubbard and spin operators (5) and (6) . It is convenient to 
collect all the terms which correspond to the same energy level (as e.g. E4 = E2, 
E5 = E3, En = E10, E14 = E12, E15 = E13) and the same number of particles 
N. Using (7) we can split H]j into 10 terms, corresponding to 10 different dimer 
energy levels (see (3)) and belonging to different subspaces of the total number 
of particles N. We obtain 

10 



Hd=Y, H d ) (8) 



1=1 

where 



H 



(i) 

D 



E2P2 + E4P4 
t x 



>?(1 



[ ) + n a 2 (l 



2^ 



(9) 
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H 



(2) 
D 



E 3 P 3 + E 5 P 5 = |[n?(l ~n a 2 -f) + n a 2 (l - n? - ^)] 



(10) 



= P 6 P 6 = - J( 2 ) • 5 2 z - ^] + ^ • ft" + ft" • 5 2 +) (11) 



H^=E 7 P r 



H 



(5) 



#8^8 = 

([/- * ) 



-^^[d+d 2 +d+d 1 ], 
— + [ AC -c\\[Si-S2- 



(12) 



<-l'<-2- 
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(t/-v) 



1 



1 + 



U — V — \ P 
2C ) + C 
U-V- J( 2 ) \ t 



2C 



+ 2C 



[d+d 2 + d+d{\ 
Ka-r^-^+n^l-n?-^)] 



1 + 



U-V + J (2) ) 
2C 



(13) 



P 



(U-V) 



+ 



4 

(U-V) 



1 



[/ - V - \ t 



2C 

U-V- 



2C 



C 



2C 



[dfd 2 + d%di] 



[nl(l-n a 2 -i-)+n b 2 (l-nl-i)} 



' _ U-V + JW) 
2C 



2 
(14) 



P^ = PioPo + E n Pn = 2J^[Si z ■ S 2 * + ^] 



(15) 



Pj? = E 12 P 12 + P 14 P 14 = (t ± g + ± J(2)) [n'n» + rfrfl 
ft + P + j(i)+j(2)) 

5 2jV2,a + K.aK'U 



(16) 



P 



(9) = E 13 P 13 + £ 15 P 15 = H + " + 4 JW + J(2) ) [nX + ngng] 

(-*+[/ + J(l) + J(2)) 



(17) 



5 



Hp ^ — Ei e Pi e — 



{U + JW+ J( 2 )) 




(18) 



2 



and n- 



n 



a . 6 



+ n! 




1,2), di (2) = 



a i(2), 1^1(2), T = c i(2),j. c i(2).t, i = 1 if i = 2 and i = 2 if i = 1. The partial 
Hamiltonians ((9), (10)), '((11)-(15)), ((16), (17)) and (18) belong to the sub- 
spaces of TV = 1,2,3 and 4, respectively. The expressions (9)-(18) are exact 
and when sum them up (see (8)) we obtain again the dimcr Hamiltonian in the 
form given by the expression (2), as it should be. The decomposition (8) of 
the dimer Hamiltonian (2) into 10 different parts (9)-(18) according to dimer 
energy levels possesses several, important advantages. First, it explicitely visu- 
alizes the important intrinsic two-site interactions, deeply hidden in the dimer 
Hamiltonian (2). Due to the fact that the formulae (9)-(18) are exact the in- 
formation about the competitive interactions within the model for a dimer is 
complete. Second, all of them are ascribed to each dimer energy level. It, how- 
ever, means that such interactions can be thermodynamically active only in the 
case when the corresponding level will be occupied by electrons. Third, we can 
see that the same types of interactions (but with different coupling constants) 
belong to quite different energy levels. It also means that they do not need to 
be thermodynamically active at the same time (it depends on the occupation 
of the particular levels) and that the resulting properties of the system depend 
on their competition. Fourth, the formulae (9)-(18) visualize the important fact 
that with the increase of the averaged number of electrons n =< N > the system 
will pass through different phases, depending on the result of the competition 
between different thermodynamically activated two-sites interactions. 

The most important two-site intrinsic interactions (leading to magnetism or 
superconductivity), presented in the expressions (9)-(18) can be devided into two 
classes. First of them belongs to magnetic interactions (ferromagnetic, antifer- 
romagnetic - it depends on the sign of the parametrs J^ 2 \U and V). Such 
interactions are present in the formulae (11) and (15) and describe Ising type 
interactions with coupling constants generated by and j( 2 \ The formula 
(11) contains also the transverse interaction between spins. The Heiscnberg type 
magnetic interactions can be seen in the first terms of the formulae (13) and 
(14), generated by more complex coupling constants, expressed by the model 
parameters , U, V and t. It is interesting to note that when = V = 
the first term in the formulae (13) and (14) describes ferromagnetic or antifcrro- 
magnetic interactions, similar to well-known t — J model (see e.g. Refs [65]- [75], 
[89]) because the coefficient ^L. ~ ^ for large | U | and is exactly the same as 
in the case of a real lattice. The first terms in the expressions (13) and (14) can 
also be considered as generalized t — J interactions, valid in a more general case 
of the model (2). The coupling constants at the first terms in (13) and (14) can 
be negative or positive. Thus, they can describe competitive ferromagnetic or 
antiferromagnetic Heisenberg interactions belonging to different dimer energy 
levels. In the case of the conventional Hubbard model ( = = V = 0) the 
coupling constants at the first terms in (13) and (14) are reduced to — -g- (13) 
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and -£7- (14) what means that the conventional Hubbard model for a dimcr de- 
scribes ferromagnetic and antiferromagnetic Hcisenberg interactions which com- 
pete together but they belong to different dimer energy levels. Let us note that 
the terms containing the products like dfd2 = b^a^ ^a2.j.£»2,t = °i ] c i i c ^-l c ^A 
and d^di = b^a^ai^bi^ = c% ^2i c hi c ^A describe the hopping of the Cooper 
pairs. Such terms, present in the second terms in (12), (13) and (14) are typical 
for the KPPK superconductivity models (cf. e.g. Refs [59]- [64]) with positive or 
negative coupling constants. Let us, however, note that this type of interactions 
are also present within the model (2) also in the case V = 0. The transver- 
sal products Si + S2~ (Si~ S*2 + ) are present in the formulae (11), (13) and 
(14). When using the second quantization we obtain S^S^ = — c i~i' c 2^j. c i,l c 2,T 

(Si 5*2 + = — c 2^t c i^J. C2 '-L Cl 't ) ano - these terms describe intersite Cooper pairs (cf. 
e.g. Ref [33] and papers cited therein). The application of the resonating va- 
lence bond approach (cf. Refs [28], [29], [52]) allows also to treat the terms like 
(Si ■ S~2 — Wl 4 " 2 ) in (13) and (14) when introducing the pairing operator / 2j i = 
^75( c 2,i c i.T — c 2,T c i,l)- For example, when we restrict ourselves to the subspace 

n\ = n% = 1 we obtain (S\ ■ S2 — j) = — f^ifiA an d sucn terms lead to super- 
conductivity (cf. Ref. [29]). It is also interesting to note that even in the case 
when = j( 2 ) = U = V = (see formulae (9)-(18)) the main competitive 
interactions (ferromagnetic, antiferromagnetic and superconducting) are always 
present when only the hopping parameter does not vanish (t 7^ 0). 

3 Conclusions 

Let us consider the case of a real lattice described by the Hamiltonian (1). We 
assume that the number of lattice points i(j) is equal to N (even number) and 
we decompose the lattice into a set of M = -y dimers described by the dimer 
index I, a (J, /?) where I(J) = 1,2, ..,M and a(/3) = 1,2. The Hamiltonian (1) 
can thus be replaced by the equivalent form 

H = J2 Hd.I + J2 tl.a;J. a c} a a Cj_ a , rT + J2 tl >a -j i pcf a a Cj,p t(T 

+ 2 S Jl.a;J,u n I,a^ n J,a,<J + \ ^ I ,a\ J,0 Ul ,a,<J n J,f3 ,<7 
IjtJ,a,a I^J,a^0,cr 
1 7 (2) 1 ^ j(2) 

1 2 2-~i ,J I a- J a n I, a ,°' n J- a , — <* ' 2 i—i I a- J 8 n I ,&,c n J,l3. — a 

- J2 V I.a;J,uc} a ,cj a ,Cj >a ,l c J,a,-\ ~ J] Vl >a ;J,pCj a ,cj a iCj^^Cj^J 

(19) 

where Hjjj is the dimer Hamiltonian given by the expression (2) where the 
lower dimer index I in the operators should be introduced (as e.g. Ci ;(7 — > Cj ; i iCT , 
etc.). The dimer Hamiltonian Hd.i in (19) can, however, be diagonalized and 
replaced by the expression (8) where again the dimer index / in the operators 
appearing in the expressions (9)-(18) should be introduced. The Hamiltonian 
(19), equivalent to (1), describes now "free" dimers (first term in (19)) and 
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their interactions (the next 8 terms in (19)). It is then evident that the Hamil- 
tonian (19) contains explicitely all the competitive two-site interactions present 
in the dimer Hamiltonian alone (see the formulae (8) and (9)-(18)) but supple- 
mented by the dimer interactions, represented by the second and further terms 
in (19). The main difference between Hamiltonians (1) and (19), both describ- 
ing the same physics, lies in the fact that the mentioned competitive magnetic 
and superconducting interactions are hidden in the Hamiltonian (1) whereas 
in the Hamiltonian (19) they appear now in a direct way. It, however, means 
that to find thermodynamical properties of the model we should simultanously 
introduce order parameters in four competitive channels (ferromagnetic, anti- 
ferromagnetic and superconducting (intrasite and intersite pairing)) where a 
special attention should be paid to the nontrivial fact that all competitive inter- 
actions are additionally scattered between different energy levels and therefore 
their activities have to be correlated with the occupation of these levels. 
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